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Abstract. We present here a simple model describing coexistence of solid and 
vapour phases. The two phases are separated by an interface. We show that 
when the concentration of supersaturated vapour reaches the dew-point, the 
droplet of solid is created spontaneously on the interface, adding to it a mono- 
layer of a "visible" size. 
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1. Introduction: Condensation phenomenon in the Ising model 

The phenomenon of droplet condensation in the framework of the Ising model was 
first described in the papers [DS1| . |DS2j . It deals with the following situation. 
Suppose we are looking at the Ising spins cr t = ±1 at low temperature 
occupying a d-dimensional box Tfj of the linear size 2N with periodic boundary 
conditions. If we impose the canonical ensemble restriction, fixing the total mean 
magnetization, 



to be equal to the spontaneous magnetization, m* (/3) > 0, then the typical config- 
uration that we see will look as a configuration of the (+)-phase. That means that 
the spins are taking mainly the values +1, while the values —1 are seen rarely, and 
the droplets of minuses in the box T$ are at most of the size of K (d) In N. We 
want now to put more —1 particles into the box T^, and we want to see how the 
above droplet picture would evolve. That means, we want to look at the model 
with a different canonical constraint: 




M N = m* (P)\T$\-b N 
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6 at > 0. It turns out that if 6/v-s are small, nothing is changed in the above picture; 
namely, if 

^ N n at 

2 ► as N — > oo, 

\ r rd I d+T 

rjvl 

then in the corresponding canonical ensemble all the droplets are still microscopic, 
not exceeding K (d) In TV in linear size. On the other hand, once 

hmini — > (J, 

N^oc \rpd I 

I n\ 

the situation becomes very different: among many ( — )-droplets there is one of the 
linear size of the order of {b^) 1 ^ d > TVW , while all the rest of the droplets arc 

d 

still at most logarithmic. Therefore 6jv ~ I^nI d+1 can ^ e canc d the condensation 
threshold , or dew-point . The behavior of the system at the threshold scale, i.e. 

d 

for b N = c\Tfj\ d+1 (1 + ojv(1)), is considered in the 2D case in [BCKj . Sharp 
description of the transition inside the threshold is considered in [HIKj . 

The above condensation picture suffers from one (largely esthetic) defect: 
both below and immediately above the condensation threshold the droplets are 
"too small to be visible" , i.e. they are of the size sublinear with respect to the 
system size. This defect was to some degree bypassed in [BSS . It is argued there 
on heuristic level, that in the low-temperature 3D Ising model in the regime when 
bpf is already of the volume order, i.e. bjy ~ vN 3 , the sequence of condensations 
happens, with "visible" results. In such regime one expects to find in the box Tjy 
a droplet T of (— )-phase, of linear size of the order of TV, having the approximate 
shape of the Wulff crystal, which crystal at low temperatures has 6 flat facets. One 
expects furthermore that the surface V itself has 6 flat facets, at least for some 
values of fejv- However, when one further increases the "supersaturation parameter" 
b n , by an increment of the order of TV 2 , one expects to observe the condensation 
of extra (— )-particles on one of the flat facets of V (randomly chosen), forming a 
monolayer m of thickness of one lattice spacing, and of linear size to be cN, with 
c > c cr it = c cr it (/?) , with c cr itN being smaller than the size of the facet. As b^ 
increases further, the monolayer m grows, until all the facet is covered by it. So one 
expects to see here the condensation of the supersaturated gas of (— )-particles into 
a monolayer of linear size ~ c cr uN, which is "visible". (Indeed, such monolayers 
were observed in the experiments of condensation of the Pb.) The rigorous results 
obtained in |BSSj are much more modest: the model studied there is the Solid-on- 
Solid model, and even in such simplified setting the evidence of appearance of the 
monolayer m of linear size is indirect. 

The purpose of the present paper is to consider another 3D lattice model, 
where one can completely control the picture and prove the above behavior to 
happen. Namely, we consider a system of ideal particles in the phase transition 
regime, and we put these phases - the vapour phase and the solid phase - into 
coexistence by applying the canonical constraint, i.e. by fixing the total number of 
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particles. We study the interface T, separating them, and we show that when we 
increase the total number of particles, the surface T changes in the way described 
above. More precisely, we show that for some values of concentration the surface 
r is essentially flat, but when the concentration increases up to the dew-point, a 
monolayer m of a size at least c cr itN appears on T, with N being the linear size 
of our system. 



2. Informal description of the main result 

In this section we describe our results informally. We will use the language of the 
Ising model, though below we treat rigorously a simpler model of the interface 
between two ideal particles phases. Ising model language makes the description 
easier; moreover, we believe that our picture holds for the Ising spins as well. 

Suppose we are looking at the Ising spins at = ±1 at low temperature /3~ 1 
in a 3D box Bm of the linear sizes RN x RN x 2N. The parameter R should be 
chosen sufficiently large in order to be compatible with the geometry of monolayer 
creation as described below. We impose (+)-boundary conditions in the upper 
half-space (z > 0), and (— )-boundary conditions in the lower half-space (z < 0). 
These (±)-boundary conditions force an interface T between the (+) and the (— ) 
phases in Vn, and the main result of the paper |D1| is a claim that the interface T is 
rigid. It means that at any location, with probability going to 1 as the temperature 
— > 0, the interface T coincides with the plane z = 0. If we impose the canonical 
ensemble restriction, fixing the total mean magnetization Mn to be zero, then the 
properties of V stay the same. 

We will now put more —1 particles into Vn) that is, we fix Mn to be 

M N = -b N = -5N 2 , 

and we will describe the evolution of the surface T as the parameter 5 > grows. 
The macroscopic image of this evolution is depicted on Figure 1. 

0. < s < s 1 

Nothing is changed in the above picture - namely, the interface T stays rigid. 
It is essentially flat at z = 0; the local fluctuations of T are rare and do not exceed 
A" In TV in linear size. 

1. S 1 < 6 < S 2 

The monolayer mi appears on T. This is a random outgrowth on T, of height 
one. Inside nti the height of T is typically z = 1, while outside it we have typically 
z = 0. 

For 8 close to S 1 the shape of mi is the Wulff shape, given by the Wulff 
construction, with the surface tension function f 2D (n) , n £ S 1 , given by 

r(n) = ^r 3 °(m) (2.1) 



m=(0,0,l) 



Here t 3D (to) , to G § 2 is the surface tension function of the 3D Ising model, the 
derivatives in (|2.1|) are taken at the point (0,0,1) € § 2 along all the tangents 
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n G S 1 to the sphere § 2 . The "radius" of rrti is of the order of N, i.e. it equals to 
T\ (5) N, and as 6 \ 6 1 we have r\ (5) \ r cr > 0. In particular, we never see a 
monolayer m of radius smaller than r cr N. 

As we explain below r cr should scale like i? 2 / 3 . In particular, it is possible 
to choose R in such a fashion that R > 2r cr or, in other words, for values of R 
sufficiently large the critical droplet fits into B^. 

As 5 increases, the monolayer trti grows in size, and at a certain moment 
S = 8 15 it touches the faces of the box Bm After that moment the shape of mi is 
different from the Wulff shape. Namely, it is the Wulff plaquette (see [SchSj ). made 
from four segments on the four sides of the RN x RN square, connected together 
by the four quarters of the Wulff shape of radius f\ (5) N. We have evidently 
f\ (S 15 ) = R/2. As S / 5 2 , the radius f\ (5) decreases to some value f\ (<5 2 ) N, 
with r 1 (5 2 ) > 0. 

II. S 2 < 6 < S 2 - 5 

The second monolayer rri2 is formed on the top of mi. Asymptotically it is 
of Wulff shape with the radius r 2 (5) TV, with r 2 (5) \ (S 2 ) as S \ S 2 , with 
r 2 (& 2 ) > 0- The first monolayer mi has a shape of Wulff plaquette with radius 
f\ (5) , which satisfies 

fi (5) = r 2 (5) . 

A somewhat curious relation is: 

r t (^ 2 ) i s strictly bigger than f\ (5 2 ) . 

In other words, the Wulff-plaquette-shaped monolayer mi undergoes a jump in its 
size and shape as the supersaturation parameter 5 crosses the value 8 2 . In fact, 
the monolayer mi shrinks in size: the radius fi (S) increases as S grows past S 2 . 

II. 5 5 2 h < S<5 3 

At the value S = 5 2 5 the growing monolayer m 2 meets the shrinking mono- 
layer mi, i.e. r 2 (<5 2 ' 5 ) = f\ (S 2 ' 5 ) = R/2. Past the value S 2 5 the two mono- 
layers rri2 C mi are in fact asymptotically equal, both having the shape of the 
Wulff plaquette with the same radius fi (5) = f 2 (5) , decreasing to the value 
f\ (5 3 ) = f 2 (i5 3 ) as 8 increases up to 5 s . 

III. S 3 < S < J 4 

The third monolayer rri3 is formed, of the asymptotic radius (5) A, with 
r 3 (S) \ r+ (6 3 ) as S \ 6 3 , with r+ (6 3 ) > 0. The radii of two bottom Wulff 
plaqucttcs fx (5) = r 2 (S) = r$ (5) decrease to the value (<5 3 ) as S decreases 
down to <5 3 , with r% (<5 3 ) > fj (<5 3 ) , so the two Wulff plaquettes mx,m2 shrink, 
jumping to a smaller area, as 8 passes the threshold value <5 3 . 

A complete investigation of the restricted Wulff variational problem (see (|7.5[) 
below) and, accordingly, a rigorous treatment of the interface repulsion phenom- 
enon which shows up on the microscopic level in all the regimes from II. 5 on is 
relegated to a forthcoming paper [IS] . For the rest of the paper we shall focus on 
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Figure 1. Creation and evolution of macroscopic monolayers on T as S grows 

the regimes 0, I and II in the context of a simplified model which we proceed to 
introduce. 

3. Our model 

We consider the following lattice model of two-phase coexistence. The 3D box 

B N = A N x {-N - 1/2, -N+ 1/2,..., AT - 1/2, N+ 1/2} 

is filled with two kinds of particles: v-particles (vapour phase) and s-particles (solid 
phase). Here Ajv is a two-dimensional RN x RN box; 

A N = {0,l,...,RN-l} 2 , 

and R is a constant, which we shall set later on to be big enough, in order to make 
our picture reacher. We have \Bjy\ = 2R 2 N 3 . Vapour i;-particles are occupying the 
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upper part of Bn, while solid s-particles - the lower part. Some sites of the box 
Bn can be empty. In our model the two phases are separated by an interface T, 
which is supposed to be an SOS-type surface; it is uniquely defined by a function 

h r :A° N ^{-N,-N + l,...,N}, 

where A° N is the interior of A^v- We assume that the interface T is pinned at zero 
height on the boundary dAjy, that is hr = on 8An. 

Such a surface V splits Bn into two parts; let us denote by Vat (r) and Sn (r) 
the upper and the lower halves. The set of configurations of our model consists 
thus from a surface T plus a choice of two subsets, a v C Vjv CH and a s C Sn (r) ; 
we have a vapour particle at a point x € Bn iff a; € cr v , and similarly for solid 
particles. 

The partition function Zn (/?) of our model is now given by 

Zn (/?) = 

e 3 q>{- ) 9|r|-(o|o- t ,| + 6|Viv(r)\o-t,| + c|o-.|+d|5j V (r)\tr,|)}. (3.1) 

(r>„ )0 - a ) 

Here |T| is the surface area of T, \a v \ is the number of vapour particles, while 
a,b,c,d are four chemical potentials. We want the two phases to be in the equilib- 
rium, so we suppose that 

e- a + e- b = e- c + e- d = e~ f , 

where the last equality is our definition of the free energy /. Accordingly, let us 
define microscopic occupation probabilities in vapour and solid states as 

p v = e f - a and p s = e f ~ c . 

To mimic the fact that the density of the solid state has to be higher, we impose 
the relation p v < p s . 

We will study our model under the condition that the total number of parti- 
cles is fixed, and in the leading order of N it is 2pR 2 N 3 , with p between the values 
p° and p s . Of course, flat interface at level zero should correspond to the choice 



More generally, given p = po + A, one expects to find T to be located approximately 
at the height £N above zero level, where £ satisfies 

e -(p s - P v ) = A. 

The above reasoning suggests that in our model the formation of macroscopic 
monolayers over flat interface should happen in canonical ensemble with total 
number of particles being fixed at 

2p R 2 N 3 + SN 2 = a N 3 + 5N 2 (3.2) 

with varying S. 
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We will denote by P the ("grand canonical") probability distribution on 
triples {T,a v ,a s } corresponding to the above partition function. Our main in- 
terest in this paper is the study of the conditional distribution of the random 
surface T, under condition that the total number of particles 



E = lav 



is fixed, i.e. the distribution 



- E 



a N 3 + SN' 



To study this conditional distribution we rely on Bayes' rule, 



(V E = a N 3 + SN 2 ^ 



P (e = a N 3 + 5N 2 


r) p(r) 




rAp(F) 



comes from volume order 



The control over the conditional probabilities P ^» rj 
local limit theorems for independent Bernoulli variables, whereas a-priori probabil- 
ities P (r) are derived from representation of T in terms of a gas of non-interacting 
contours. 

In the sequel C\ , c%, . . . are positive constants which appear in various in- 
equalities and whose values are fixed in such a way that the corresponding bounds 
hold true. 



4. Volume order limit theorems 

r ) is easy, since we are dealing 



The study of probabilities Pr (E = a N 3 + SN 2 

with independent variables. Indeed, let Bpj = SVU Vat be the decomposition of Bpj 
induced by T. Then, the P ^» | F^j -conditional distribution of the overall number 
of particles is 

s = E e? + E S> 

with iid Bernoulli (p s ) random variables £f , and iid Bernoulli random variables 

^3 

Let a(T) be the signed volume under the interface F, 

a(T) = J J h T {x,y)dxdy, (4.1) 

where we set hr to be equal to hr(i) in the unit box i + [1/2, 1/2] 2 . Clearly, 
\S N \ = R 2 N 3 + a(T) and [Vzvj = R 2 N 3 - a(T). Accordingly, 

E(E r) = a Q N 3 + a(T)p sv , 

where p sv = p s — p v . Introducing the variances D s = p s (\ — p s ), D v = p°(l — p v ) 
and D = D s + D v , we infer from the Local Limit Theorem (LLT) behavior: For 
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every K fixed there exist two positive constants C\ and C2, such that 



ci < 



a N 3 + 5N 2 



■■ exp 



<C 2 , 



(4.2) 



^ttD\B n \ 

uniformly in N, \S\ < K and T, provided |a(r)| < KN 2 . 



5. Surface weights 

We now want to describe the a-priori probability distribution P (T) . It is convenient 
and natural to express it via the weights {w (T)}, so that 

P ( r)4p r( r) = ^_ > (5.i) 

where we shall use an additional symbol Pr in order to stress that the corresponding 
probabilities are computed in the contour model we are going to introduce now. 

For our purposes it is necessary to introduce a contour parameterization of 
the set of all surfaces T. Contours will live on the bonds of the dual (two di- 
mensional) box = {1/2, 3/2, . . . , RN — 3/2} , and they are defined as follows: 
Given an interface T and, accordingly, the height function hr which, by definition, 
is identically zero outside A^, define the following semi- infinite subset T of M 3 , 

?= |J ((x,y,k) + d) , 

where C = [—1/2, 1/2] 3 is the unit cube. The above union is over all (x,y) e 1? 
and k e 1/2 + Z. 

Consider now the level sets of T, i.e. the sets 

H k = H k (rj = {(x,y)eR 2 :(x iy ,k)€?}, k = -N, -N + 1, . . . , N. 

We define contours as the connected components of sets dHt- The length | — ^ | of a 
contour is defined in an obvious way. Since, by construction all contours are closed 
polygons composed of the nearest neighbour bonds of A^, the notions of interiour 
int(7) and exteriour ext(7) of a contour 7 arc well defined. A contour 7 is called a 
©-contour (O-contour), if the values of the function hr at the immediate exterior 
of 7 arc smaller (bigger) than those at the immediate interiour of 7. 

Alternatively, let us orient the bonds of each contours 7 C dHk in such a 
way that when we traverse 7 the set Hk remains to the right. Then ©-contours 
are those which are clockwise oriented with respect to their interiour, whereas 
©-contours are counter-clockwise oriented with respect to their interiour. 

Let us say that two oriented contours 7 and 7' are compatible, 7^7', if 

1. Either int(7) n int(7') = or int(7) C in.t(V) or int(V) C int(7). 
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2. Whenever 7 and 7' share a bond b, b has the same orientation in both 7 and 



A family T = {7.;} of oriented contours is called consistent, if contours of 
r are pair-wise compatible. It is clear that the interfaces T are in one-to-one 
correspondence with consistent families of oriented contours. The height function 
hr could be reconstructed from a consistent family r = {7} in the following way: 
For every contour 7 the sign of 7, which we denote as sign(7), could be read from 
it orientation. Then, 



where xa is the indicator function of A. 

We are finally ready to specify the weights w(T) which appear in (|5 . ip : Let 
r = {7} be a consistent family of oriented (signed) contours, Then, 



By definition the weight of the flat interface w(ro) = 1. 
6. Estimates in the contour ensemble 

In order to make the contour model (|5.ip , (|5.2|) tractable one should, evidently, 
make certain assumptions on the largeness of (3, e.g. e* 3 should be certainly larger 
than the connective constant of self-avoiding random walks on 1? [MS] . In fact, it 
would be possible to push for optimal results in terms of the range of /? along the 
lines of recent developments in the Ornstein-Zcrnike theory [II ICIVH ICIV2] . How- 
ever, in order to facilitate the exposition and in order to focus on the phenomenon 
of monolayer creation per se, we shall just conveniently assume that (3 is so large 
that one or another form of cluster expansion goes through, see eg. }D2j . Due to 
the (±-contour) symmetry of the model the corresponding techniques would be 
quite similar to those developed in the context of the 2D low temperature Ising 
model in [DKSj . Consequently, instead of stating conditions on f3 explicitly we 
shall just assume that f3 > /3o, where (3q is so large that all the claims formulated 
below are true. 

In the sequel we shall employ the following notation: C for clusters of non- 
compatible contours and 3>g(C) for the corresponding cluster weights which shows 
up in the cluster expansion representation of partition functions. 

Peierls estimate on appearance of"/. Given a contour 7 and a consistent family 

k 

of contours T, let us say that 76 T, if 7 appears in T exactly k times. Then, 



h-y(x,y) 




■yer 




(5.2) 




(6.1) 
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Indeed, every V satisfying 7 6 V can be decomposed as T = T' U 7 U • ■ ■ U 7. 
Therefore, 

£ r , u;(r')e- fe/3 lTl 



Pr her 



< 



where the sums are over all consistent families which are compatible with 7, but 
do not contain it. 

Fluctuations of a(T) and absence of intermediate contours. The following rough 
a-priori statement is a consequence of (|6.1[) : There exist positive v such that for 
every bo > fixed, 



Pr(|a(r)| > bN 2 ) < c 3 e 



(6.2) 



uniformly in b > bo and in N large enough. 

In view of (|4.2p (computed with respect to the flat interface Tq with a(To) 



0) the bound (|6.2p implies that the canonical distribution P f < 
is concentrated on T with 



a(r) < TV 2 max 



v 2 D 2 R 4 



Jk 



E = a iV 3 + (57V 2 
(6.3) 



Now let the interface T be given by a consistent collection of contours, and assume 
that 7 ~ T. Of course a(r U 7) = a(T) + a(j). Let us assume that the surface T 
satisfies the estimate (16.31). Then 



(ru 7 



S = a N J 



SN 2 



< 



P (S = a iV 3 + 5A^ 2 


ru 7 ) 


Pr (r U 7) 


P = a iV 3 + <5iV 2 


r ) 


Pr(r) 



< C4 exp < C5 



N7)l 



/3|7l ?■ < c 4 cxp <^ c 6 



/3|7l 



where we have successively relied on Bayes' rule, ()4.2|) and on the isopcrimctric 
inequality 

It follows that for every K there exists e = e((3) > such that intermediate 
contours 7 with 

-\ogN < I7I < eiV (6.4) 
e 

are, uniformly in |5| < if, improbable under the conditional distribution 

S = a N 3 + SN 2 ^ 



In the sequel we shall frequently ignore intermediate contours, as if they do not 
contribute at all to the distribution (|5.ip . To avoid confusion, we shall use Pr for 
the restricted contour ensemble, which is defined exactly as in (|5.1|) . except that 
the intermediate contours 7 satisfying (|6.4|) are suppressed. 
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7. The surface tension and the Wulff shape 

Since we are anticipating formation of a monolayer droplet on the interface, we 
are going to need the surface tension function in order to study such a droplet and 
to determine its shape. It is defined in the following way: Let A be an oriented 
site self avoiding path on the dual lattice Z^. An oriented contour 7 is said to be 
compatible with A; 7 ~ A, if A n int(7) = and if whenever A and 7 share a bond 
b, the orientation of b is the same in both A and 7. Accordingly, if C is a cluster of 
(incompatible) contours, then C ~ A if 7 ~ A for every 7 G C. 

In the sequel 0* = (1/2, 1/2) denotes the origin of ZjJ. Let x G Z^. Set, 

Tp(x) = J2 ™p\-P\M - £<M C > 

A:0*^x [ C/A 

where the sum is over all oriented self-avoiding paths from 0* to x. 

Let n € S 1 be a unit vector, and n 1 - G S 1 is orthogonal to it. The surface 
tension Tp in direction n is defined as 

T P {n) = - lim \\ogTp{[Ln^\). 

Consider the Wulff variational problem, which is a question of finding the 
minimum wp (S) of the functional, 

wp (S) = min W (A) . 

{A:Aroa(A)=S} 

Here 

W (A) = J Tfs (n s ) ds, 

n s being the unit normal to A at the point A (s) , and the minimum is taken over all 
closed self-avoiding loops A, enclosing the area S. Of course, wp (S) — V~Swp (1) . 
Let us denote by Wp the Wulff shape, which is the minimizing loop with area 
5=1. 

As in [DKSj it could be shown that if (3 is sufficiently large, then rp is well 
defined and strictly positive. Furthermore, the boundary of the optimal loop Wp 
is locally analytic and has uniformly positive and bounded curvature. 

One can now apply to the present setting the machinery and the results of 
[DKS] , [DST] , [DS2] , [Sets] and (iSch) . They allow us to study the probabilities of 
the events 

Pr(A 6 ) = Pr{F:a(r) = b} , (7.1) 

where we consider here the probability distribution (|5.ip . 

As it follows from local limit results in the restricted phase |DKSj without 
intermediate contours (|6.4p . for all values of b, the probability Pr (Ab) is bounded 
above by 

Pr (At) < c 7 exp |-c 8 -^ A n\ . (7.2) 
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In particular, for the values of b <C A 3 / 2 the main contribution to Pr (At,) comes 
from small contours; |7| < (T 1 log N. In other words, for such values of b, condi- 



tional distribution Pr I 



is concentrated on the interfaces T which are essen- 
tially flat: all contours 7 of a typical surface T are less than e _1 log TV in length, 
while their density goes to zero as (3 — > 00. 

On the other hand, for values of b ^> A 3 / 2 long contours contribute, and the 
probabilities Pr (Ab) satisfy 

logPr(A 6 ) = -Vbwp (1) (1 + o N (1)) , (7.3) 

provided, of course, that the scaled Wulff shape ^Jb/N 2 Wp fits into the square 
[0, R} 2 . Under these two restrictions on b the analysis of |DKSj implies that the 

conditional distribution Pr (■ A^j is concentrated on the interfaces T which arc 

"occupying two consecutive levels". Namely, the set {7^} of contours, comprising 
r, contains exactly one large contour, 70, of diameter ~ s/b, while the rest of them 
have their lengths not exceeding e _1 In A. The contour 70 is of ©-type, so for the 
majority of points inside 70 the value of the height function hp is 1, while outside 
70 it is mainly zero. Finally, the contour 70 has 

Asymptotic shape: The contour 70 is of size ~ y/b, and it follows very close 
the curve VbWp. Namely, the latter can be shifted in such a way that the 
Hausdorff distance 

Ph (70, VbW ) < #b. (7.4) 



Of course, all the claims above should be understood to hold only on the set 
of typical configurations, i.e. on the sets of (conditional) probabilities going to 1 
as A — > 00. 

In the present paper we also we need to consider such values of b ~ 2R 2 N 2 , 
when the scaled Wulff shape ^Jb/N 2 Wp does not fit into the square [0, R} 2 . This 
situation was partially treated in the paper [SchSj , and the technique of that paper 
provides us with the following information about the typical behavior of T under 

the distribution Pr (■ A^j for the remaining values of b. 

Namely, instead of the Wulff variational problem we have to consider the 
following restricted Wulff variational problem, which is a problem of finding the 
minimum 

w r / t (S)= i min W5* t (A;Ai,...,A fc )=W(Ai) + ... + W(A fc ), (7.5) 
where 

• the curves Ai, A& are closed piecewise smooth loops inside the unit square 

Qx\ 

• the loops Ai are nested: Int (Afe) C Int (Afe_i) C ... C Int (Ai) ; 

• Area (A fe ) + Arca(A fc _i) + ... + Arca(Ai) = S. 

The parameter k is not fixed; we have to minimize over k as well. For the area 
parameter S small enough, the minimum in (|7.5|l is attained at k = 1, while Ai is 
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the scaled Wulff shape, s/SWp. In other words, in this regime Wp st (5) — uup (5) . 

Let Si be the maximal value, for which the inclusion y/SWp C Qi is possible. In 
the range Si < S < 1 the solution to (|7.5|) is given by k = 1, while the loop Ai is 
the corresponding Wulff plaquette, described above. In the range 1 < S < 25i the 
solution has the value k = 2, the curve Ai is the Wulff plaquette, while the curve 
A2 C Ai is the Wulff shape; they are uniquely defined by the two conditions: 

1. Area (A 2 )+ Area (Ai) = 5, 

2. the curved parts of Ai are translations of the corresponding quarters of A2. 

In the range 25i < S < 2 we have k = 2, while the loops A2 = Ai are identical 
Wulff plaquettes. 

The relation (|7.3[) is generalized to 

logPr (A b ) = -RNw^ (jpjp) (1 + o N (1)) . (7.6) 

The function w r p st (5) is evidently increasing in S. For S small it behaves as 

dy/S. In the vicinity of the point S = 1 it behaves as d' \J 5 — 1 for 5 > 1, and 
as d"\fl — 5 for 5 < 1. Otherwise it is a smooth function of 5, < 5 < 2. The 
two singularities we just pointed out, are responsible for the interesting geometric 
behavior of our model, which has been described informally in Section 2, and will 
be explicitcly formulated in the next Section. Namely, each one is responsible for 
the appearance of the corresponding droplet. 

Accordingly, once the Wulff shape \JbjN 2 Wp docs not fit into the square 
[0, R] 2 , while b < c (/?) R 2 N 2 ( where the constant c (/?) — > 1 as (3 — > 00) the condi- 
tional distribution Pr A^j is concentrated on the interfaces T which again arc 

occupying two consecutive levels. The set {7^} of contours, comprising T, contains 
one large contour, 70, this time of diameter ~ R, which in some places is going 
very close to the boundary of our box. The rest of contours have their lengths not 
exceeding e~ J IniV. The contour 70 is of ©-type, and for the majority of points 
inside 70 the value of the height function hr is 1, while outside 70 it is mainly zero. 
Finally, the contour 70 has asymptotic shape of the Wulff plaquette, in the same 
sense as in (|7.4[) . 

In the remaining range R 2 N 2 < b < 2R 2 N 2 the set {7^} of contours, com- 
prising r, contains exactly two large contours, 70 and 71, with 71 C Int (70) , both 
of the ©-type. The interface T is, naturally, occupying three consecutive levels: it 
is (typically) at the height 2 inside 71 , at height 1 between 70 and 71 , and at height 
outside 70. Note that for b close to R 2 N 2 the contour 71 is free to move inside 
70, so its location is random (as is also the case in the regime of the unique large 
contour, when the scaled Wulff shape yJbjN 2 Wp fits into the square [0, R] 2 )- The 
contour 70, on the other hand, is (nearly) touching all four sides of the boundary 
of our box, so it is relatively less free to fluctuate. 

In the complementary regime, when b is close to 2R 2 N 2 , the two contours 70 
and 71 have the same size in the leading order (which is linear in N), while the 
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Hausdorff distance between them is only ~ iV 1 / 2 ; it is created as a result of the 
entropic repulsion between them. In particular, in the limit as N — ► oo, and under 
the scaling, the two contours coincide, going in asymptotic shape to the same 
Wulff plaquette. The study of this case needs the technique, additional to that 
contained in |DKS| , [DSlj , |DS2j , [SchSj and [ISch] , since the case of two repelling 
large contours was not considered there. The case of the values b above 2R 2 N 2 
is even more involved, since there we have to deal with several large mutually 
repelling contours. We will return to it in a separate publication, see [IS] . 



8. Main result 

We are ready now to describe the monolayers creation in our model: Let us fix 
p s > p v (and hence p sv and D), and let (3 be sufficiently large. Let us also fix R 
large enough, so that the rescaled Wulff shape of area 



lD 2 w 2 (1) 



R 4 / 3 



fits into the R x R square. 



Theorem 8.1. Let T be a typical interface drawn from the conditional distribution 
P ( • I E = a N 3 + SN 2 ) . Define 



S 1 = -yD 2 w 2 p sv R 4/3 . (8.1) 

For values of S satisfying < 5 < S 1 , the interface T is essentially flat: all 
contours ofT have lengths bounded above by e~ 1 \ogN. 

There exists S 2 > 6 , such that for 6 1 < 5 < S 2 the interface T has one 
monolayer. Precisely, T contains exactly one large contour 70 of approxi- 
mately Wulff shape (or Wulff plaquette shape), such that 

«(7o)>^^ 2 - (8-2) 

The rest of contours of T are small; their lengths are bounded above by 
e" 1 log TV. 

Similarly, there exists a value 5r, such that for 5 2 < S < 5r the interface T 
has two monolayers, and contains exactly two large contours, 70 and 71 C 
Int(7o). The bigger one, 70, has the shape of the Wulff plaquette, while the 
smaller one has the Wulff shape. Again, 0(71) > N 2 . 



9. Proof of the main result 

Let us fix S and consider the surface distribution P ( • | E = ao + SN 2 ) . Since 
we can ignore intermediate contours (|6.4[) and since we already know how the 
typical surfaces looks like in the constraint ensembles Pr ( • | ^4^) i h would be 
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enough to study conditional probabilities P ( A b | S = ao + SN 2 ). Namely, for ev- 
ery 5 we need to know the range of the typical values of the "volume" observ- 
able b. To do this we will compare the probabilities P(yl&, S = ao + <5iV 2 ) ~ 
P ( £ = ao + (57V 2 1 Ab) Pr (A;,) for various values of b, in order to find the domi- 
nant one. 

There are three regimes to be worked out: Fix r\ € (0,1/2) and eg small 
enough. 

Case 1. b < N 1+r >. By ([42|) and fT2l) . 



( S 2 ( b 2 

< iV 3/2 P (E = a + SN 2 | A) Pr (Aj,) (9.1) 



2 



S 2 Sv sv b b 2 

Z = a + 5N 2 \A b )Pr(A b ) < c 12 exp - — -.N + -f^ - c 8 -, A N 



Case 2. iV 1+ '' < b < c 9 N 2 . By g^) and 

2DR 2 ^' ' NR 2 D " a N 2 

(9.2) 

Obviously, once c g is chosen to be sufficiently small, the right hand side of (|9.2p 
is negligible with respect to the lower bound on left-hand side of (|9.ip (computed 
at 6< iV 1+r '). 

Case 3. b = pN 2 with p > c 9 . By {715) and, once again, by volume order local limit 
result 

— p SV p) 2 _ npr rst f P\ 



exp { dr 2 N RNw e vm ) - 0{N) 

< P (S = a + SN 2 | A b ) Pr (A 6 ) (9.3) 

^- i ^~-«-?'(£)+H 

Therefore, in order to figure out the dominant contribution between (19. ip 
and (|9.3p , we have to find the global minimum of the function 

^W(£) (M) 

on the interval p £ [0, 2i? 2 ]. This minimization problem needs just the elementary 

calculus, see e.g. |BCKj . For small values of p our function reduces to ^j^r — l~ 
wp (1) y/p. After the following change of variables: 



VP 

A = — — and n = k(S) 



^3/2 



2DR 2 wp (l)v / P 7; 
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A) 2 + VA. 




3/2 



(9.5) 



One easily sees that 

• If k < k c , then the global minimizer is 0. 

• If k = k c then there are exactly two global minimizers; and A c = 2/3. 

• If k > k c , then the global minimizer A m is the maximal solution of 

4nV\(l - A) , 
which, in particular, satisfies A m > 2/3. 

A similar analysis applies in the vicinity of the singularity of the function 
w r p st (jp) at ~ 1. Since the function w r p st (S) is monotone, and has the deriv- 
ative equal to +00 at S = 1, the point of the global minimum of (|9.4p . which is a 
monotone function of <5, never belongs to some neighborhood of the point = 1. 
Therefore at some 5 = 5 2 it jumps from some value p_ < R 2 to p + > R 2 . 

The proof of Theorem 1 is, thereby, completed. 

10. Conclusions 

In this paper we have described a model of the interface between the vapour and 
liquid phases, evolving as the total number of particles increases. We have shown 
that the evolution of the interface goes via the spontaneous formation on it of 
one monolayer of the size of the system. We believe that the same result can be 
proven for the 3D Ising model with the same boundary conditions, i.e. periodic 
in two horizontal directions and ± in the vertical one. It will be very interesting 
to establish the phenomenon of the monolayer formation in the 3D Ising model 
with (+)-boundary conditions, when the monolayer attaches itself to a facet of the 
Wulff-like (random) crystal. This problem, however, seems to be quite difficult, 
since one needs to control the rounded part of the crystal. This rounded part is 
probably behaving as a massless Gaussian random surface (compare with [K]). 
and this alone indicates enough the complexity of the problem. 
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